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Abstract 

Models of supersymmetric grand unified theories based on the SU(5)gut x U(./V)h 
gauge group (N = 2, 3) have a symmetry that guarantees light Higgs doublets and 
the absence of dimension-five proton-decay operators. We analysed the proton decay 
induced by the gauge-boson exchange in these models. Upper bounds of proton lifetime 
are obtained; r(p — > ir°e + ) ^ 6.0 x 10 33 yrs in the SU(5)gutxU(2)h model and r(p —* 
7r°e + )<5.3 x 10 35 yrs in the SU(5)gut x U(3)h model. Various uncertainties in the 
predictions are also discussed. 



1 Introduction 



The supersymmetric (SUSY) grand unified theories (GUTs) are among one of the most 
promising candidates of physics beyond the standard model (SM); this is because of their 
theoretical beauty and is also because the gauge-coupling unification is supported by precision 
experiments. The GUTs generically predict a proton decay through a gauge-boson exchange 1 . 
The lifetime of the proton, however, varies very much from model to model; it is proportional 
to the fourth power of the gauge-boson mass. Thus, the proton-decay experiments are not 
only able to provide strong support for the GUTs but also to select some models out of 
several candidates. 

For more than twenty years, many attempts have been made at constructing models of 
SUSY GUTs. There are two strong hints to find realistic models. First, the two Higgs 
doublets are light, whereas their mass term is not forbidden by any of the gauge symmetries 
of the minimal supersymmetric standard model (MSSM). Second, the rate of proton decay 
through dimension- five operators [2] is smaller than naturally expected [3 2 . Thus, the two 
phenomena require two small parameters; it would then be conventional wisdom to consider 
that there might be symmetries behind the small parameters. Moreover, only one symmetry 
is sufficient to explain the two phenomena because of the following reason. If a symmetry 
forbids the mass term of the two Higgs doublets, 



then the dimension-five proton-decay operators are also forbidden by the same symmetry, 



and vice versa. Here, it is implicitly assumed that the quarks and leptons have the same 
charge under the symmetry when they are in the same SU(5) multiplet, and that the Yukawa 
couplings of quarks and leptons are allowed by the symmetry. 

There are three classes of models that have such a symmetry 3 . One uses the SU(5)i 
xSU(5) 2 gauge group [TQ1 EJ where there is an unbroken Z N symmetry pHH2]. The 
second class consists of models based on the SU(5)gut><U(./V)h (N = 2,3) gauge group 

1 There is a class of models of SUSY GUTs where the gauge-boson exchange does not induce the proton 
decay 0. 

2 The experimental lower bound of the proton lifetime through those operators is r(p — ► K + v) ^ 6.7 x 10 32 
yrs (90% C.L.) 0. 

3 There are other types of models that are not constructed as field theoretical models in four-dimensional 
space-time 00,0 and 00. Reference is a string-theoretical realization of 03 EUE2J an d Refs. 00 
are higher-dimensional extension of ^31- The proton decay is analysed in for the model and in [To] 

for the model 0. 
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[TBI H3 E] > where there is an unbroken discrete R symmetry [T3] . The last class of models 
can be constructed so that the unified gauge group is a simple group [TH]. The symmetry 
discussed in the previous paragraph, however, cannot remain unbroken; it is broken in such 
a way that the dimension-five operators are not completely forbidden. The proton decay 
through the gauge-boson exchange is discussed in [Tj5J |20] for these models. 

In this article, we analyse the proton decay for the second class of models. The dimension- 
five operators are completely forbidden in these models, and the proton decay is induced by 
the gauge-boson exchange. Our analysis is based on models in [T3|, which use the SU(5)gut x 
U(-/V) H gauge group (N = 2, 3). Reference [2TJ obtained an estimate of the proton lifetime in 
the SU(5)gut x U(3)h model, adopting a number of ansatze to make the analysis simple. The 
estimate was 4 r(p — > 7r°e + ) ~ (0.7-3) xlO 34 yrs, and hence there is an intriguing possibility 
that the proton decay is observed in the next generation of water Cerenkov detectors. This 
article presents a full analysis: both the SU(5)gut x U(2)h model and the SU(5)gut x U(3)h 
one are analysed without any ansatze. Three parameters of the models are fixed by three 
gauge coupling constants of the MSSM, and the remaining parameters are left undetermined. 
The range of these parameters is restricted when the models are required to be in the calcu- 
lable regime. As a result, the range of the gauge-boson mass is restricted. Thus, we obtain 
the range 5 of the lifetime rather than an estimate of it. 

The organization of this article is as follows. First, we briefly review both the SU(5)gutX 
U(A^)h models (N = 2, 3) in section[21 The range of the GUT-gauge-boson mass is determined 
both for the SU(5)gutxU(2)h model and for the SU(5)gutxU(3)h model in sections El and 
EJ respectively. In particular, it is shown that the range of the mass is bounded from above, 
which leads to an upper bound of the lifetime of the proton in each model. The upper bounds 6 
are, in general, predictions that can be confirmed by experiments. Since the SUSY-particle 
spectrum affects the MSSM gauge coupling constants through threshold corrections, the 
upper bound of the lifetime depends on the spectrum. Therefore, the upper bound is shown 
as a function of SUSY-breaking parameters in sectionElfor both models. Various uncertainties 
in our predictions are also discussed. Section [HI gives a brief summary of the results obtained 
in this article, and compares the results with predictions of other models. We find that 
r(p -> vr°e + ) < 6.0 x 10 33 yrs in the SU(5) G utxU(2) h model, and r(p -> 7r°e + ) < 5.3 x 10 35 
yrs in the SU(5)gut x U(3)h model; here, we exploit the uncertainties in the value of the 
QCD coupling constant (by ±2cr) and in the threshold corrections from the SUSY particles; 
other uncertainties, which cannot be estimated, are not included in these figures. 

4 The following numerical value is based on a re-analysis in |22| . 

5 This procedure is the one adopted in [2Sj, where the minimal SU(5) SUSY GUT model was analysed. 
6 Only a lower bound is obtained, e.g. in the minimal SU(5) model |2.'ij . 
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2 Brief Review of Models 



2.1 SU(5) gut xU(2)h Model 

Let us first explain a model based on the product gauge group SU(5)gutxU(2)h- Quarks 
and leptons are singlets of the U(2) H gauge group and form three families of 5*+10 of the 
SU(5) gut- Fields introduced to break the SU(5)gut symmetry are given as follows: X 01 ^ 
(a, (3 = 1,2), which transforms as (l,adj.=3 + 1) under the SU(5)gutxU(2)h gauge group, 
and Q a i (i = 1,...,5) + Q% and Q\ (i = 1,...,5) + Q 6 a , which transform as (5*+l,2) and 
(5+1,2*). The ordinary Higgs fields H l (5) and Hi(5*) are not introduced; the fields Q$, 
Q & a are eventually identified with the two Higgs doublets. The SU(5)gut index is denoted 
by i, and the U(2)h index by a or (3. The chiral superfield X a g is also written as X c (t c ) a /3 
(c = 0, 1, 2, 3), where t a (a = 1, 2, 3) are Pauli matrices of the SU(2) H gauge group 7 and t = 
l 2x2 /2, where U(2) H — SU(2) h xU(1)h- The (mod 4)-R charge assignment of all the fields in 
this model is summarized in Table ^ This symmetry forbids both the enormous mass term 8 
W = QsQa and the dangerous dimension-five proton-decay operators W — 10 • 10 • 10 • 5*. 
The most generic superpotential under the R symmetry is given by 

W = V2X 2 uQ i a X a (t a ) a f3 Q l3 i + V2\' m Q 6 a X a (t a ) a l3 Q l3 6 
+V2\ m Q\X (t r^ + V2\> m Q 6 a X%t r^ 6 

-V2X m v 2 X a a (3) 
+c 10 10 ill2 10 i3i4 (gQ) i5 6 + c 5 .(QQ) 6 4 • W ij ■ 5* + • • • , 

where the parameter v is taken to be of the order of the GUT scale; A2H, A' 2H , Am and A' 1H are 
dimensionless coupling constants; C10 and 05. have dimensions of (mass) -1 . Ellipses stand 
for neutrino-mass terms and other non-renormalizable terms. The fields and Q l a in the 
bifundamental representations acquire vacuum expectation values (VEVs), (Q a i) = t><5" and 
(Q l a) = v ^ l a j because of the first three lines in (jSJ). Thus, the gauge group SU(5)gut x U(2) h 
is broken down to that of the SM. The first terms in both the first and second lines in (j3J) 
provide mass terms for the unwanted particles. Q$ and Q & a are identified with the Higgs 
doublets in this model, and no Higgs triplets appear in the spectrum. As a result, no particle 
other than the MSSM fields, not even a gauge singlet of the MSSM, remains in the low-energy 

7 The normalization condition tr(i a £fc) = <S a b/2 is understood. Note that the normalization of the following 
to is determined in such a way that it also satisfies tr(toio) = 1/2. 

8 The mass term of the order of the weak scale for the Higgs doublets can be obtained through the 
Giudice-Masiero mechanism 24 . 
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spectrum. This fact not only guarantees that the gauge coupling unification of the MSSM 
is maintained, but also that the vacuum is isolated. The R symmetry is not broken at the 
GUT scale, and the /x term, dimension-four and dimension-five proton-decay operators are 
forbidden by this unbroken symmetry. 

Fine structure constants of the MSSM are given at tree level by 

1 1 1 



a 3 a c ctcuT 



(4) 



1111 

— = — = + , (5) 

a 2 OIL «GUT «2H 

and 

1 3 1 3 

— = = — — + (6) 

a x a Y a G UT «ih 

where ckgut, «2H and are fine structure constants of SU(5)gut, SU(2) h and U(1)h, 
respectively. Thus, the approximate unification of a 3 , a 2 and a\ is maintained when a 2 n and 
am are sufficiently large. 

Although it is true that the gauge coupling unification is no longer a generic prediction 
of this model, nevertheless, it need not be a mere coincidence. Here, the gauge coupling 
unification is a consequence of the fact that a 2 n and am are relatively strongly coupled 
when compared with acur- 

2.2 SU(5) GUT xU(3) H Model 

The other model is based on an SU(5)gutxU(3)h gauge group, where U(3)h — SU(3)hX 
U(1)h- Under the SU(5)gutxU(3)h gauge group, the particle content of this model is Q k a (5+ 
1,3*), Q a k {5* + 1,3) and X° /3 (l,adj=3<g>3*), with k = 1, . . . , 5, 6; a, [3 = 1, 2, 3. In addition, 
there are the ordinary three families of quarks and leptons (5*+10,l) and Higgs multiplets 
H l +Hi (5+5*, 1). An R symmetry forbids both (1) and (2). The R charges of the fields are 
summarized in Table El 

The most generic superpotential under the R symmetry is given JH] by 

W = \/2X w Q i a X a {t a ) a pQ p i + V2X' m Q 6 a X a (t a )°j 3 Q l3 6 
+V2X lu Q\X°(t r^ l + V2X' m Q 6 a X°(t r^ 6 
-V2X m v 2 X a a (7) 

+y 10 10 ■ 10 • H + y 5 *5* • 10 •# + ••• , 
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where t a (a = 1, 2, 8) are Gell-Mann matrices, t = 1.3x3/^, Vio and t/ 5 * are Yukawa cou- 
pling constants of the quarks and leptons, and A3H, A' 3H , Aih, A' 1H , b! and h are dimensionless 
coupling constants. The first three lines of (J7J) lead to the desirable VEV of the form (Q a j) 
= vS^ and (Q l a ) = v5 l a . Thus, the SU(5)gut><U(3)h gauge group is broken down to that 
of the SM. The mass terms of the coloured Higgs multiplets arise from the fourth line in (J7J) 
in the GUT-symmetry-breaking vacuum. No unwanted particle remains massless. 

The fine structure constants of the SU(3)hxU(1)h groups must be larger than that of the 
SU(5)qut- This is because the gauge coupling constants of the MSSM are given by 



and 



1 


1 


1 


+ 


1 


a 3 


ac 


CtGUT 


"3H 


1 


1 


1 








a L 


Q-GUT 


1 
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3 

= _5_ 
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2 
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ai 


ay 


a-GUT 


am 



where «3H and am are fine structure constants of SU(3)h and U(1)h, respectively. Thus, the 
approximate unification of 0-3, a 2 and a\ is maintained when a 3 n and «ih are sufficiently 
large. 



2.3 Rough Estimate of Matching Scale 

Figure ^ shows the renormalization-group evolution of the three gauge coupling constants of 
the MSSM. The tree-level matching equations (jlj)-© and (P)- (fTU|) suggest that the matching 
scale is below the energy scale M 2 _3 in Fig. [T]in the SU(5)gutxU(2)h model, and is between 
the two energy scales M 2 _3 and Mi_ 2 in the SU(5)gutxU(3)h model. Here, M 2 _3 is the 
energy scale where coupling constants of the SU(2)x and the SU(3)c are equal, and Mi_ 2 
the one where coupling constants of the SU(2) L and the U(l)y are equal. In particular, the 
matching scale is lower than the scale Mx_ 2 , i.e. the conventional definition of the unification 
scale, in both models. Thus, the decay rate of the proton is higher than the conventional 
estimate, which uses Mi_ 2 ~2x 10 16 GeV as the GUT-gauge-boson mass. 
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3 Gauge-Boson Mass in the SU(5) gut xU(2)h Model 



Let us now proceed from the discussion at tree level to the next-to-leading-order analysis 
in order to draw more precise predictions. To find the proton-decay rate, it is necessary to 
determine the mass of the gauge boson, rather than the matching scale. The GUT-gauge- 
boson mass enters the threshold corrections to the gauge coupling constants at the 1-loop 
level, and it can hence be discussed directly. 

The analysis in this article follows the procedure described in [23] . First, the three gauge 
coupling constants of the MSSM are given in terms of gauge coupling constants and other 
parameters of GUT models. We include 1-loop threshold corrections from the GUT-scale 
spectra to the matching equations. Then, we constrain the parameters of GUT models by 
the three matching equations: three parameters are determined, and other parameters are 
left undetermined. The free parameters, however, cannot be completely free when we require 
that the GUT models be in a calculable regime, i.e. when a perturbation analysis is valid. 
We determine the calculable region in the space of the free parameters and, as a result, the 
ranges of the GUT-gauge-boson masses are obtained for the models. 

In the minimal SU(5) SUSY-GUT model, for example, there are four parameters. The 
three gauge coupling constants of the MSSM determine the coloured-Higgs mass, and put 
two independent constraints between the other three parameters. The three parameters 
are the unified gauge coupling constant, the GUT-gauge-boson mass and a coefficient of 
the cubic coupling of the SU(5)cuT-adj. chiral multiplet in the superpotential [23]. Thus, 
the matching equations cannot determine the GUT-gauge-boson mass directly. The cubic- 
coupling coefficient is chosen as the free parameter, while the GUT-gauge-boson mass and 
the unified gauge coupling constant are solved in terms of the free parameter and the MSSM 
gauge coupling constants. The free parameter, however, cannot be too large; otherwise it 
would immediately make itself extremely large in the renormalizat ion-group evolution toward 
the ultraviolet (UV). Thus, it is bounded from above, and its upper bound leads to the lower 
bound of the GUT-gauge-boson mass of the minimal SU(5) model [25] . 

The SU(5)gutxU(./V)h model (N = 2 (or 3)) has five (or six) parameters in the three 
matching equations of gauge coupling constants, as we see later. Thus, two (or three) pa- 
rameters are left undetermined. The space of two (or three) free parameters is restricted 
by requiring the perturbation analysis to be valid, just as in the analysis of the minimal 
SU(5) model. As a result, the range of the GUT-gauge-boson mass is obtained. The crucial 
difference between the three models is that only the lower bound of the mass is obtained in 
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the minimal SU(5) model, while the upper bound is obtained both in the SU(5)gut x U(2) h 
model and in the SU(5)gutxU(3)h model 9 , as shown in the following. The SU(5)gut><U(2)h 
model is analysed in this section, and the result of the SU(5)gut x U(3)h model is described 
in section HJ 

3.1 Parameters of the Model 

The MSSM gauge coupling constants are given in terms of parameters of the SU(5)gutX 
U(2)h model at the 1-loop level as 

a 3 a GV T 2ir \MJ 2?r \MJ ' 

J-0.) = _J_ w + J_ w + zi + A ln (^) + _L ln (^K), (12) 

«2 «GUT «2H 27T \M J 2-k \ M J 2n \M 3G J 

-0*) = — M + J-W + ^lnf^ + ^lnf^V (13) 



Oil «GUT CKih 27T \M/ 27T \ M 

where M and /i are renormalization points of the GUT model and MSSM, respectively. The 
renormalization point M is chosen to be above the GUT scale and fi to be just below the GUT 
scale. The right-hand sides consist of the tree-level contributions (the first and second terms) 
in Eqs. (0])-©, 1-loop renormalization and threshold corrections (the remaining terms). 
Gauge coupling constants are considered to be defined in the DR scheme, and hence the 
step-function approximation is valid in the 1-loop threshold corrections [23]. Various mass 
parameters of the model enter the equations through the threshold corrections; Mq is the 
GUT-gauge-boson mass, M 3V and M 3C are masses of the SU(2)^-adj. vector multiplet and 
chiral multiplet, respectively. These mass parameters are given in terms of the parameters 
of the Lagrangian (at tree level), as shown in Table El 

There are five parameters of the GUT model in the above three equations: Mq, M 3 v/M 3 c, 
1/o>2e, 1/qih and I/ckqut- Three of them are solved in terms of the other two parameters and 
of the three MSSM gauge coupling constants. The other two parameters are left undetermined 
for the moment. We take Mq and M 3 y/M 3 c as the two independent free parameters. Then 
three others, namely aQUT(M G ), a 2 n(M G ) and aw(M G ), are determined through Eqs. (jTTjl - 
(JT3J) by setting \i = M = M G . Another parameter of the model, «2 H = (A 2 h) 2 /(47t), is also 
expressed in terms of a 2 n{M G ), aQUT(M G ) and M 3V /M 3G : 

1 (m o) j * wj^y. (i4, 



a 2H ' \otw{M G ) + a GVT (M G 



3 The lower bound also exists in this model. 
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3.2 Parameter Region of the Model 



Let us now determine the parameter region in the parameter space spanned by Mq and 
M 3 v/M 3 c- We require that the perturbation analysis be valid; it is necessary that all the 
coupling constants in the model be finite in the renormalization-group evolution toward the 
UV, at least within the range of spectrum of the model. To be more explicit, the coupling 
constants aza(M), a^iM), a^^M), axa(M), a^ H (M) and a^ H (M) are required to be finite 
in the renormalization-group evolution, at least while the renormalization point M is below 
the heaviest particle of the model. We use this necessary condition to determine the parameter 
region. Here and hereafter, we adopt the following notation: 

«2H - ^ , «2H - An , «1H - , «1H - An ■ W 

First, we use the 1-loop renormalization-group equation 10 to determine the parameter 
region. Renormalization-group equations of this model are listed in appendix El The result 
is shown in the left panel of Fig. El The parameter region is given by the shaded region 
in the Mq-(M 3 v/M 3 q) plane. The analysis is based on the value of ol 3 {h) calculated from 
a MS,(5) (Mz) = 0.1212, i.e. the value larger than the central value by 2a. The reason for this 
choice is explained shortly. 

The result in the left panel is understood intuitively as follows. First, the parameter Mq 
is bounded from above (from the right of the panel). It is quite a natural consequence, since 
it is consistent with the rough estimate of the matching scale in subsection 12.31 Secondly, 
the parameter region is bounded also from below. It is also a natural consequence because 
of the following reason. The beta function of the superpotential coupling ct^n in Eq. (jlUj) 
implies that this coupling constant immediately becomes large unless its contribution to the 
beta function is cancelled by those from gauge interactions. Thus, the parameter space with 
«2H *C a^, which is almost equivalent to M 3V <C M 3C , is excluded. 

We adopt the value a^ s, ^(Mz) = 0.1212 for the value of the QCD coupling constant, 
rather than the usual central value a^ s '^(Mz) = 0.1172, because this allows M2-3 to be 
larger. In turn, this allows the excluded region to be smaller, and hence the upper bound for 
the GUT-gauge-boson mass becomes more conservative with this choice. 

10 The beta function of a^^M) depends also on a^nfM), a^ H (M) and a^ H (M). Thus, the beta function 
cannot be calculated without the values of those coupling constants. Their values, however, are not deter- 
mined through the 1-loop matching equations (f X 1|) — (|13|) . Therefore, we set them in the beta function as 0, 
so that a2n(M) becomes large as slowly as possible in the evolution to the UV. This makes the excluded 
parameter space smaller and makes our analysis more conservative. 
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We further include 2-loop effects in the beta functions of the gauge coupling constants 11 . 
The renormalization-group equations for the gauge coupling constants are listed in appendix 
1X1 The 2-loop effects also become important at generic points of the parameter space, because 
the 1-loop beta functions of the gauge coupling constants are accidentally small everywhere 
on the parameter space 12 . The initial values of the 2-loop renormalization-group evolution, 
i.e. values at the matching scale M = Mq, are not determined for o^hj a m an d ce m from 
the matching equations (fTTj) - (fT3|) . Thus, we set their values as 

«2H = a m, «m = «ih = «ih, (16) 
when the renormalization point M is at Mq. Although we should have varied also these 
values as free parameters, we believe that the result of our analysis is not affected very much 
by changing these values; the reason is explained in appendix [BJ The right panel of Fig. El is 
the parameter region determined in this analysis. 

The right panel of Fig. El shows that the parameter space with M 3 y 3> M 3 c, i.e. « 2 h 3> 
a 2H , * s further excluded, and the only surviving parameter region is around the line of M%y ~ 
M3C, i.e. «2H — «2H- fr 1S clear, as shown below, why this region and only this region 
survives. Let us neglect, for the moment, the renormalization effects from the SU(5)gut 
gauge interaction; the SU(5)gut gauge coupling constant is smaller than those of the SU(2) H 
and the U(1)h interactions. Then, one can see that the 2-loop part of the beta functions of 
a 2 n and am are proportional to (a 2 H ~ a m) an d (&m ~ a m)- Thus, the renormalization 
effects from « 2 h and am are completely cancelled 13 by a^ and a m just in that region. 

The cancellation in the 2-loop beta functions is due to the (A/ - = 2)-SUSY structure in the 
GUT-symmetry-breaking sector [HI El EI] ; the beta functions of gauge coupling constants 
vanish at two loops and higher in perturbative expansion in the M = 2 SUSY gauge theories 
j2H|. Therefore, the remaining region at the 2-loop level survives even if higher- loop effects 
are included in the beta functions. 

The renormalization-group evolution is determined by the 1-loop beta functions on the 
(J\f = 2)-SUSY line « 2 h — « 2 h (when the SU(5)gut interaction is neglected). Therefore, 
we consider that the point in the parameter space indicated by an arrow in the right panel 
of Fig. El gives a conservative upper bound of Mq- We also consider that the upper bound 
so obtained is a good approximation of the maximum value of Mq, although the parameter 

n Note that the beta functions are scheme-independent up to two loops for gauge coupling constants, while 
only up to one loop for coupling constants in the superpotential. 

12 At 1-loop order, the beta function of a2n is given as a result of cancellation between vector loop 3Tg = 6 
and chiral loop — = —8. 

13 Here, we assume that Eqs. (tTB|) are also satisfied. 
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region becomes thinner and thinner as Mq increases; see appendix |B]for a detailed discussion. 
Theoretical uncertainties in this determination of the conservative upper bound of Mq are 
discussed in subsection 13.31 A related discussion is also found in appendix |Bj 

Now that we know that the upper bound is obtained on the [M = 2)-SUSY line, it is 
possible to obtain the upper bound of M G without numerical analysis. Indeed, the following 
two facts make the analysis very simple; Mq is essentially the only free parameter on the 
line, and the 1-loop renormalization-group evolution is a good approximation there. 

The gauge coupling constant « 2 h is given at Mq by 

(17) 



a 2H (M, 




G) 

through the matching equations (JTTj) and ()12j) . where a threshold correction proportional 
to ln(M 3 y/M 3 c) is neglected owing to the Af = 2 SUSY. Here, M 2 _ 3 is defined so that 
a 3 (M 2 _ 3 ) = a 2 (M 2 _ 3 ). The gauge coupling constant a 2 n so determined should not be too 
large because 

1 ~ 1 Afcf^W ai 



a m {M w ) a m {M G ) 2tt V M G 

It follows only from the inequality in (fTHj) 14 that (27r)/a 2H (M G ) £ 3.7; note that M 3V /M G 
can be expressed in terms of am and «gut- Thus, the upper bound of Mq is given by 

M G ~ e ~ ia ^ M a) M 2 _ z < e"^M 2 _3 ~ 0.40 x M 2 _ 3 . (19) 

3.3 Uncertainties in the Upper Bound of the Gauge-Boson Mass 

Here, we estimate uncertainties in our prediction of the upper bound of the GUT-gauge- 
boson mass. Uncertainties arising from our analysis of the GUT model are discussed in this 
subsection, while uncertainties arising from low-energy physics are discussed in section |S] 

First, we focus on the effects from the SU(5)gut interaction. They have been neglected 15 
in the discussion of the previous subsection, but they do contribute to the 2-loop beta func- 
tions; in addition, the higher-loop contributions from a 2 n and ai^ no longer cancel because 
the SU(5)gut interaction does not preserve N = 2 SUSY. Thus, the renormalization-group 
evolution is changed and the determination of the upper bound is affected. The SU(5)gut 
interaction contributes to the beta function of am in Eq. ()39|) by less than 10% of the 1-loop 

14 "3.7" is almost independent of the SUSY threshold corrections to the MSSM gauge coupling constants. 
15 They are not neglected in the numerical analysis in Fig. 
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contribution 16 . Thus, the value of «2h(^g) f° r the upper-bound value of M G is not changed 
by 10% (see Eq. |T8]) ). As a result, the upper bound of Mq is not modified by a factor of 
more than e (2T)/(4a 2H (M G ))x(±io%) ^ io± - 04 . 

Second, the perturbative expansion would not converge when the 't Hooft coupling 
2a 2 H/(4vr) exceeds unity. It is impossible to extract any definite statement on the renormali- 
zation-group evolution when the perturbative expansion is not valid. However, most of the 
renormalizat ion-group evolution is in the perturbative regime, i.e. (47r/(2a2H)) ~ 1 ; since we 
know that (47r / (2« 2 h)) (Mg) — 3.7 for the upper-bound value of M G . Thus, we consider that 
the perturbation analysis in the previous subsection is fairly reliable. 

Third, non-perturbative contributions are also expected in the beta functions, and they 
might not be neglected since the gauge coupling constants are relatively large in this model. 
They 17 are expected to be of the form [27] 



oo 



^^(^^^(^r, (20) 



where c n 's are numerical factors of the order of unity. Each contribution comes from n- 
instantons. Here, we neglected perturbative and non-perturbative contributions through 
wave-function renormalizations of hypermultiplets. This is because hypermultiplets of M = 
2 SUSY gauge theories are protected from any radiative corrections [2S]. We see that the 
non-perturbative effects given above are not significant when the renormalization point is 
around the GUT scale, since 



-4 2 ^ 



(02h) e «2H ~5x 1(T 5 < 1. (21) 

So far, the analysis is based on a renormalizable theory. However, the Yukawa couplings of 
quarks and charged leptons are given by non-renormalizable operators in Q. Another non- 
renormalizable operator is also necessary to account for the fact that the Yukawa coupling 
constants of strange quark and muon are not unified in a simple manner. Those operators, 
however, affect the gauge coupling constants through renormalization group only at higher- 
loop levels. Moreover, they are not relevant to the renormalization-group flow (say, in the 
sense of Wilsonian renormalization group) except near the cut-off scale M*. These are the 
reasons why we neglected the effects of those operators. 
16 See appendix IbI for more details. 

17 We thank Tohru Eguchi for bringing this issue to our attention. 
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There may be, however, a non-renormalizable operator 



W = 2tr\\± + c^\W°W a \, (22) 

which directly modifies the matching equations of the gauge coupling constants at tree level. 
Exactly the same analysis as in subsections 13. II and 13.21 tells us that the upper bound of M G 
given in Eq. (19) is modified 18 into 

M G < 0.40 x io-°- 82xc x M 2 _ 3 , (23) 

as long as c>(— 3.7/(67r) + 14/9 x (e 9 /a G UT))/l-2 ~ -0.16 + 1.3 x (e 9 /a G uT); here, e g = 
((93 -gi)/gi)(M^ 2 ) ~ -(0.03-0.01). 



4 Gauge-Boson Mass in the SU(5) GUT xU(3) H Model 
4.1 Parameter Region of the Model 

The same analysis as in section El is performed for the SU(5)gut x U(3)h model. The 1-loop 
matching equations in this model, which are quite similar to Eqs. ()lip -([13 p . are found in [2*T] . 
The particle spectrum around the GUT scale, which comes into the threshold corrections, is 
summarized in Table 0] 

There are six parameters in the matching equations: M G , M SV /M 8C , I/ogut, l/ a 3H, and 
l/ a iH, just as in the SU(5)gut x U(2) h model, and M Hc M H -/M G . Three of them are fixed 
through the matching equations, and the other three are left undetermined. We take M G , 
M 8V /M 8C and M Hc M H -/M G as the three free parameters. The space of these parameter is 
restricted by requiring that all the coupling constants am(M), a 3 ^ , ain(M) and be 
finite in the renormalization-group evolution toward the UV, at least within the range of the 
spectrum. 

The parameter region is shown in Fig. El only the \J (Mh c Mjj 5 / M G ) = 10 a3 cross section 
is described, and hence the region is described in the M G -(M 8 v / M 8G ) plane. This analysis is 
based on the value of a^(n) that is calculated from a^ s '^ B ' (M z ) = 0.1132, i.e. a value smaller 
than the central value by 2a. This is because it makes our analysis more conservative. The 
parameter region obtained by the 1-loop renormalization group is shown as the shaded area in 

18 Contributions to Eqs. (|llfl H3|) are, for example, c x I6ir x 0.024 ~ 1.2 x c at the point in the parameter 
space indicated by an arrow in Fig. [21 
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Fig. |21 The region is bounded from the right and from the left, which is again consistent with 
the rough estimate of the matching scale given in subsection 12.31 The region is also bounded 
from below just for the same reason as in subsection 13.21 The parameter region distant 
from the (Af = 2)-SUSY line 19 M 8 y ~ M 8 c is excluded when 2-loop effects are included 20 in 
the beta functions of the gauge coupling constants; the 2-loop contributions have significant 
effects compared with 1-loop effects, because the 1-loop beta function of accidentally 
vanishes. The 1-loop renormalization-group evolution is reliable on the (A/" = 2)-SUSY line, 
the thick line in Fig. [31 and hence the points indicated by (A) and (B) give the upper and the 



lower bound of Mq, respectively, in the y (Mh c Mh e /M g ) = 10 3 cross section. The upper 
and the lower bound of Mq of the model are the maximum and the minimum value that Mq 
takes at (A) and (B), respectively, when the remaining parameter (Mh c Mh-/M g ) is varied. 
Since it is evident from the figure that the lower bound of Mq leads to too fast a proton 
decay, we focus in the following only on the upper bound of Mq. 

4.2 Uncertainties in the Upper Bound of the Gauge-Boson Mass 

In this subsection, we estimate uncertainties in the prediction of the upper bound of the GUT- 
gauge-boson mass obtained in the previous subsection. The uncertainties that originate from 
low-energy physics, however, are discussed in section 

First, we discuss the effects of the interactions that violate the M = 2 SUSY. The 
SU(5)gut gauge interaction and the cubic couplings in the fourth line of (J7J) are the sources 
of the violation of the M = 2 SUSY. Those interactions change the 1-loop-exact evolution 
of M = 2 SUSY gauge theories. The change in the upper bound of Mq comes 21 from the 
change in the evolution of am, because the upper bound was determined by the running of 
am in the absence of [M = 2)-SUSY breaking. The beta function of am is changed at most 
by a few per cent 22 , which leads to the change of the upper bound of Mq by a factor of at 



19 Thc TV = 2 SUSY is enhanced in the GUT-symmetry-breaking sector when 

5m ^ A m (~ A' 1H ), g m ~ A 3H (~ A 3H ), (24) 

are satisfied and h, h , « GU t are neglected [T71I51I2T]. 

20 We set the initial values (M = Mq) of coupling constants that are not determined by the matching 
equation as follows: , , \ \' \ \' 

_ h = h', al H = a$ n , ^ H =^„=a 1H . (25) 

This choice makes the renormahzation-group evolution the most stable. 

21 We confirmed that the change in the evolution of is not so significant as to make the finiteness of 

a more severe condition than that of am- 

22 This estimate comes from the ratio between the 1-loop contribution and the SU(5)gut contribution 
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most 10 ±001 . 

Second, one can see from the matching equations [2Tj of this model that the gauge coupling 
constants a 3 u and am are not so large as to invalidate the perturbative expansion when the 
coloured-Higgs particles are moderately heavier than the GUT-gauge-boson mass; only one 
threshold correction from the coloured-Higgs particles is sufficient to keep both coupling 
constants within the perturbative regime. Non-perturbative effects are not important at all 
in such a region. 

Finally, a non-renormalizable operator that corresponds to (|22j) may also exist in this 
model. Such an operator, if it exists, contributes to the matching equations at tree level. In 
its presence, we can perform exactly the same analysis as in the previous subsection. The 
result of this analysis is presented in section |BJ 



5 Conservative Upper Bound of Proton Lifetime 

The analysis in sections EU and [3] presented the way of extracting the upper bound of the GUT- 
gauge-boson mass for both models. The lifetime of the proton through the GUT-gauge-boson 
exchange is given [23| in terms of Mq as 



n , , « /0.015GeV 3 \ /2.5\ 2 / 1 \ / M G \ 4 

r (p-> vr°e + )~ 1.0xl0 35 x — — yrs, 26 

yP } \ a n J \A R J v25a GU T(M G ) J Uo 16 GeW y ' V } 

where an is a hadron matrix element 23 calculated with lattice quenched QCD [3U] (an = 
—0.015 ±0.001 GeV 3 ) renormalized at 2.3 GeV and Ar ~ 2.5 a renormalization factor of the 
dimension-six proton-decay operators [3T]; Ar consists of a short-distance part, A^^ ~ 2.1, 
which comes from the renormalization between the GUT scale and the electroweak scale, and 
a long-distance part 24 , A^ D ^ ~ 1.2, from the renormalization between the electroweak scale 

at two loops. Note also that the ah = /i 2 /(47r) contribution has a sign opposite to that of the SU(5)gut 
contribution. 

23 The hadron matrix element «h is defined by (vac.|(uijdij)ui|j5(k)) = a^u(k.). This is related to another 
matrix element W (~ -0.15 ±0.02 GeV 2 ) through 

an 



V2U 



-W, (27) 



where W is defined by lim (7r (p)|(uii;(ifl)ui|p(k)) = Wu(k) , and / w is the pion decay constant (2/ ff = 130±5 
MeV E2|). 

24 The numerical coefficient of the formula of the lifetime adopted in [21] is different from the one in Eq. I|2()l) 
in this article. This is because the formula in |2] is based implicitly on Ar ~ 3.6 in [22], whose value is the 
effect of renormalization between the GUT scale and 1 GeV. It was therefore incorrect, in [2J, to use at the 
same time Ar ~ 3.6 renormalized at 1 GeV and the hadron matrix element in |30j renormalized at 2.3 GeV. 
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and 2.3 GeV (Ar = A R ■ A R ). We note the expression of A R [32] f° r later convenience: 

4/3 3/2 23/30 

, ( sd) = ( ac(M z ) \ "^ ( a L (M z ) \ ( a Y (M z ) \ ^ 
R \a c (M G )) [a L (M G )J \a Y (M G )) ' 1 ) 

where 6, (i = 1,2,3) are coefficients of the 1-loop beta functions of the three gauge coupling 
constants of the MSSM. The renormalization from Yukawa coupling constants is omitted 
because its effect is negligible. 

Threshold corrections from SUSY particles are of the same order as those from the par- 
ticles around the GUT scale. The 2-loop effects in the renormalization-group evolution 
between the electroweak scale and the GUT scale are also of the same order. Therefore, the 
two above effects should be taken into consideration in deriving predictions on the GUT- 
gauge-boson mass (and hence on the lifetime of the proton). This implies, in particular, that 
the predictions depend on the spectrum of SUSY particles. We present the predictions on 
the upper bound of the lifetime of the proton as a function of SUSY-breaking parameters 
of the mSUGRA boundary condition in subsection 15.11 Predictions can be obtained also 
for other SUSY-particle spectra such as that of gauge-mediated SUSY breaking (subsection 
15. 2|) . Subsection 15.31 discusses how the predictions are changed when there are vector-like 
SU(5)GUT-multiplets at a scale below the GUT scale. 

5.1 mSUGRA SUSY Threshold Corrections 

Let us first consider the SUSY-particle spectrum determined by the mSUGRA boundary 
condition. This spectrum and the MSSM gauge coupling constants in the DR scheme are 
calculated in an iterative procedure. We use the S 'OFTS US 'Yl. 7 code [33! for this purpose. 
These coupling constants are evolved up to the GUT scale through the 2-loop renormalization 
group. They are used as input in the left-hand sides of, say, Eqs. (fTTj) - (fT3*j) . to obtain a 
prediction of the upper bound of the GUT-gauge-boson mass. The universal scalar mass mo 
and the universal gaugino mass are varied, while we fix the other parameters of the 
mSUGRA boundary condition 25 ; tan/3 = 10.0, A = GeV, and the sign of the \i parameter 
is the standard one. 

The left panel of Fig. |U is a contour plot on the m -m 1 / 2 plane, describing the upper 
bound of the proton lifetime in the SU(5)gutxU(2)h model, where we set the unknown 
coefficient c of the non-renormalizable operator (22) to zero. The QCD coupling constant 
25 This is because changes in these parameters did not change the result at all as in |21j . 
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a MS,(5) (Mz) = 0.1212 is used, so that the upper bound becomes more conservative. One 
can see that this upper bound ranges over (1.4-3.2) xlO 33 yrs. Notice that the (m , m^) 
dependence of the proton lifetime arises almost only through the variation of M 2 _3 (see Eq. 
(Unj). Indeed, the contours of the upper bound of the lifetime in the left panel of Fig. |U 
behave in the same way as those of M 2 _3 in the upper-left panel of Fig. 

It is now easy to see how much the prediction is changed when we adopt the central value 
of the QCD coupling constant, af s, ^(M z ) = 0.1172. Since the choice of the QCD coupling 
constant directly changes M2-3, it severely affects the upper bound of Mq in this model. 
M 2 _3 is decreased by a factor e ^-^ l o.n72 0.1212^ anc i ^he lifetime is shortened by a factor 
e -27rxo.28 ^ q 2. We confirmed that the upper bound of the lifetime does not exceed 1 x 10 33 
yrs even when the SUSY-breaking parameters m and mi/2 are varied up to 2000 GeV if we 
adopt the central value of the QCD coupling constant. 

The hadron matrix element an in [3U], which has a statistical error an = —0.015 ± 
O.OOlGeV 3 , does not include a systematic error (e.g. an error due to the quenched approxi- 
mation). Reference (35] estimates that the systematic error is about 50%, which leads to an 
uncertainty in the lifetime of a factor of 2. 

Therefore, the conservative upper bound is roughly r^ 6 x 10 33 yrs, where we exploit the 
uncertainties in the SUSY threshold corrections, in the value of the QCD coupling constant 
and in the hadron matrix element. Thus, the prediction does not contradicts the experi- 
mental lower bound from Super-Kamiokande r(p — > 7r°e + )> 5 x 10 33 yrs (90% C.L.) |H1 at 
this moment 26 , yet the large portion of the parameter region is already excluded. Moreover, 
one can expect that the uncertainties originating from low-energy physics will be reduced 
in the future. Thus, further accumulation of data in Super-Kamiokande and the next gen- 
eration of water-Cerenkov detectors will be sure either to exclude this model without the 
non-renormalizable operator (|2*2*|) . or to detect the proton decay. 

Now, we move to consider the SU(5)gutxU(3)h model. The right panel of Fig. 0] is a 
contour plot on the mo-mi/2 plane, describing the upper bound of the proton lifetime in the 
SU(5)gutxU(3) h model. The QCD coupling constant c^ 5 )(M z ) = 0.1132 is used. The 
upper bound of the proton lifetime ranges over (1-5) xlO 35 yrs on the mSUGRA parameter 
space that is not excluded by the LEP II bound on the lightest-Higgs-boson mass. 

Let us now see how much the above prediction is changed by uncertainties related to 
the QCD. First, the following observation is important in discussing the effect from the 

26 The lifetime listed in is r{p -> 7r°e+) > 1.6 x 10 33 yrs (90% C.L.), based on a paper jSZ] published in 
1998. 
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uncertainty in the value of the QCD coupling constant. The behaviour of the contours of 
the upper bound of the lifetime, and hence of the GUT-gauge-boson mass has, in this model, 
strong correlations with that of M\_ 2 presented in the upper-right panel of Fig. We find 
an empirical relation 

M G <0.60 x Mi_ 2 . (29) 

Thus, the upper bound does not depend on the value of the QCD coupling constant very- 
much, since Mi_ 2 is not affected very much. Second, the uncertainty in the hadron form 
factor is common to both models. Therefore, the most conservative upper bound of the 
proton lifetime is roughly r ^ 10 36 yrs in this model. In particular, the proton decay might 
not be within the reach of the next generation of experiments. 

5.2 Threshold Corrections from Various SUSY-Particle Spectrum 

Gauge-mediated SUSY breaking (GMSB) is one of the highly motivated models of SUSY 
breaking. The spectrum of the SUSY particles is different from that of the mSUGRA SUSY 
breaking, and moreover, there are extra SU(5)cuT-charged particles as messengers. Thus, 
the predictions on the proton lifetime are different from those in the case of mSUGRA 
SUSY breaking. We discuss the effects of the difference in the SUSY-particle spectra in 
this subsection. A possible change of predictions due to the existence of extra particles is 
discussed in the next subsection. 

The ranges of the GUT-gauge-boson masses is different for different SUSY-particle spec- 
tra, yet the difference only arises from the difference in the two energy scales M 2 _ 3 and 
Mx_ 2 : the energy scale where the SU(2) L and SU(3)c coupling constants become the same 
and where those of U(l)y and SU(2)^ become the same, respectively. The upper bound of 
Mq is given in terms of M 2 _3 through Eq. (|T9*j) in the SU(5)gutxU(2)h model, and in terms 
of Mi_ 2 through Eq. flU in the SU(5) G utxU(3) h model. 

Figure shows how M 2 _ 3 and M x _ 2 vary over the parameter space of the GMSB. The 
parameter space is spanned by two parameters: an overall mass scale A of the SUSY breaking 
in the MSSM sector and the messenger mass M mess . We assume that the messenger sector 
consists of one pair of SU(5)gut-(5+5*) representations. Gaugino masses are given by 

*"* = S a ( 1 + °(m^)) (30) 

at the threshold M mess . We calculate the SUSY-particle spectrum, the SUSY threshold 
corrections to the MSSM gauge coupling constants and the renormalization-group evolution 
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to the messenger scale using the code (32] • We include contributions from the messenger 
particles into the beta functions in the renormalization-group evolution from the messenger 
scale to the GUT scale. M 2 _ 3 and Mi_ 2 are obtained and are shown in Fig. It is clear 
from Fig. that the ranges of M 2 _3 and Mi_ 2 are almost the same in the mSUGRA and 
GMSB parameter spaces. Therefore, we conclude that there is little effect that comes purely 
from the difference between the SUSY-particle spectra of the mSUGRA and GMSB. 

The gaugino masses satisfy the GUT relation in both the mSUGRA and GMSB spectra, 
which may be the reason why M 2 _ 3 and M x _ 2 are almost the same in the two spectra. The 
gaugino mass spectrum, however, might not satisfy the GUT relation 27 . Even in this case, 
we can obtain the upper bound of the lifetime through M 2 _3 for the SU(5)gutxU(2)h model 
and through Mi_ 2 for the SU(5)gut><U(3)h model. 

5.3 Vector-Like SU(5)GUT-Multiplet at Low Energy 

There are several motivations to consider charged particles in vector-like representations, 
whose masses are of the order of the SUSY-breaking scale or an intermediate scale. Messenger 
particles are necessary in the GMSB models, and the anomaly cancellation of the discrete R 
symmetry also requires [39J extra particles such as SU(5)gut-(5+5*). 

There are three effects on the proton lifetime in the presence of these particles. The first 
two effects come from the change in the values of the unified coupling constant ogut and the 
renormalization factor Ar of the proton-decay operators. First, the unified gauge coupling 
constant is larger in the presence of new particles, and hence the decay rate is enhanced. Then, 
the lifetime is shortened by a factor not smaller than 0.66 when a vector-like pair SU(5)gut- 
(5+5*) exists at an energy scale not lower than 1 TeV. Second, the renormalization factor 
Ar is changed by such a vector-like pair only in its short-distance part. The new expression 
for A^ D ^ is now given by 



(SD) 



4/3 t > 3/2 , ^ 23/30 

33/5 



' ac{M z ) \ 3 f a L (M z ) \ - 1 ( a Y {M z ) \ 

v a c (M) ) \ a L (M) J { a Y (M) J 

4/3 3/2 
/ / a c (M) \ 2 / a L (M) \ 2 ( a Y (M) ' 

s \a c {M G )) \a L {M G )) \a Y (M G ) 



23/30 
-38/5 

(31) 



where M is the mass scale of the vector-like pair. We find that A^^ increases from 2.1 to 2.5 

27 Gaugino masses without the GUT relation are not unnatural at all in the product-group unification 
models we discuss in this article |38| . 
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as the mass scale M decreases from the GUT scale to 1 TeV. Thus, the lifetime is shortened 
by a factor not smaller than 0.71 because of the renormalization factor. 

The third effect is due to threshold corrections from the vector-like particles. The triplets 
and doublets in the vector-like pair 5+5* are expected to have different masses, just as the 
bottom quark and tau lepton do. The triplets will be heavier than the doublets by 

M 3+3 » f aciMs+^Y ] 1 (32) 
M 2+2 * ^ a c (M G ) J ' 1 ] 

which increases from 1.0 to 2.1 as the mass scale M5+5* of a vector-like pair SU(5)gut- 
(5+5*) decreases from the GUT scale to 1 TeV. The upper bound of the proton lifetime 
in the SU(5)gutxU(2)h model becomes tighter by a factor (M 2+ 2«/M 3+ 3») ^0.48 as M 2 _ 3 
is decreased by a factor (M 2+2 * / 'M 3+3 ») 1//4 ; in the SU(5)gut x U(3)h model, instead, it is 
loosened by a factor (M 3+3 */M 2+ 2*) 2 ^ 7 ~ 1-2 as M x _ 2 is increased. 

The proton decay is, thus, enhanced by all three effects in the SU(5)gutxU(2)h model; 
the lifetime is shortened by a factor of 0.22 when SU(5)gut-(5+5*) exists at 1 TeV. The rate 
is enhanced also in the SU(5)gutxU(3)h model; the lifetime is shortened by a factor of 0.56. 



6 Conclusions and Discussion 

We analysed the proton-decay amplitude in a class of models of SUSY GUTs: SU(5)gut 
xU(]V)h models with N = 2,3. Dimension-five proton-decay operators are completely for- 
bidden, and hence the gauge-boson exchange is the process that dominates the proton decay. 
We found that the gauge-boson mass is bounded from above by 

M G < 0.40 x i(r a82xc x M 2 _ 3 (33) 

in the SU(5) G utxU(2) h model 28 and by 

M G < 0.60 x i(T a4xc x Mi_ 2 (34) 

in the SU(5)gut><U(3)h model. Here, M 2 _ 3 (Mi_ 2 ) denotes an energy scale where SU(2) l and 
SU(3)c (U(l)y and SU(2)^) gauge coupling constants are equal, respectively. In the right- 
hand sides, c are coefficients of non-renormalizable operator (|22j) in the SU(5)gut x U(2)h 
model and of the one that corresponds to (J22J) in the SU(5)gut x U(3)h model. It is quite 
28 This expression for the upper bound of Mq is valid as long as —0.16 + 1.3 x (£ s /q;gut)- 
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important to note that the upper bound was obtained in these models (for fixed c), which 
leads to the upper bound of the lifetime. Although the gauge-boson masses are bounded 
also from below in the latter model, the lower bound is of no importance. This is because it 
predicts a lifetime much shorter than the lower bound obtained so far from experiments. 

The coefficients c directly affect the gauge coupling unification, and hence they appear 
in the above formulae. One will see later that they are the largest source of uncertainties 
in the upper bound of the lifetime if c are of the order of unity. Although there may be an 
extra (broken) symmetry or any dynamics that suppress the non-renormalizable operators, 
we leave c in the formulae for generic cases. 

In section ^ we briefly mentioned two other classes of models of SUSY GUTs constructed 
in four-dimensional space-time. Let us make a brief summary on the mass of the gauge 
bosons of such models before we proceed to a discussion of the lifetime. 

Let us first discuss the gauge-boson mass in the models in (TUJ [TT] . The spectrum around 
the GUT scale consists of three ((adj.,l)° + (l,adj)°)'s and two ((3,2)" 5 / 6 + h.c.)'s of the 
MSSM gauge group, in addition to the GUT gauge boson. Parameters of the models allow 
a spectrum where the matter particles are lighter than the GUT gauge boson. Then, 1-loop 
threshold corrections from such a spectrum imply that the GUT-gauge-boson mass is heavier 
than the energy scale of approximate unification [13] ■ Therefore, no upper bound is virtually 
obtained in the models in ^UJ [nj. A lower bound might be obtained, yet no full study 
has been done so far. Non-renormalizable operators in the gauge kinetic functions affect the 
matching equations just as in our analysis. 

On the contrary, in the models in non-renormalizable operators do not affect the 
matching equations and, moreover, the mass of the GUT gauge boson is smaller than the 
energy scale where the three gauge coupling constants are approximately unified: 

M G ~ A a M unif , (35) 

where M uni f ~ M 2 _3 ~ Mi_ 2 , A a small parameter of the order of 10 _1 and a the charge of 
a field whose VEV breaks the SU(5)gut symmetry down to SU(3)cxSU(2)iXU(l)y. Thus, 
Mc is fairly small in the models. The upper bound would be obtained once a model (A and 
a, in particular) is fixed. The Super-Kamiokande experiment already puts constraints on the 
choice of A and a. The proton decays also through dimension-five operators, although these 
operators can be suppressed in some models in this class. 

Thus, the ranges of the proton lifetime of those models lie in the following order: 

r([IHl)di m .6 ~ r(SU(5) GUT x U(2) H ) < r(SU(5) GUT x U(3) H ) < r(HD|, HH). (36) 
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However, the ranges would have certain amount of overlap between one another, and hence 
it would be impossible to single out a model only from the decay rate of the proton. Detailed 
information on the branching ratio of various decay modes does not help for that purpose 
either; the decay is induced in all the above models 29 by one and the same 30 mechanism: the 
gauge-boson exchange. 

Even if no model can be singled out, one can, and one will be able to exclude some of the 
models on the basis of experimental results currently available and obtained in the future, 
respectively. We summarize, in the following, the upper bound of the proton lifetime for the 
SU(5)gut x U(2)h model and the SU(5)gut x U(3)h model. It would also be of importance 
if one finds an upper bound and a lower bound of the lifetime in models in jTHj and in 
PHD!! H21, respectively. 

Now the proton lifetime is bounded from above by 



r(^ 7 r°e + )<4.1xl0 32 ^ 3 (l(T 3 - 2c ) UUXJ ^ V l±± yrs 

^ ; ~ U0 158 GeW V J [ au J \25a GVT J \A B J y 

(37) 

in the SU(5)gut x U(2)h model. The largest uncertainty in this prediction comes from the 
value of c. Another from the systematic error in an. No estimate is available for the value 
of c. The error in «h is not studied very much, yet the lifetime is changed by a factor of 
(0.5-2) if the conservative estimate in |35j is adopted 31 . The experimental value of the QCD 
coupling constant a^ s '^(Mz) changes the prediction through the change in M 2 _3. The upper 
bound is changed by a factor of (0.15-5.9) when the coupling constant is varied by ±2a error 
determined by experiments. The threshold corrections from SUSY particles also changes the 
prediction through the change in M2-3. They change M2-3 typically from 10 15,77 GeV to 
10 15 ' 90 GeV, and hence the upper bound is changed by a factor of (0.75-2.5). Therefore, the 
theoretical upper bound exceeds the experimental lower bound (r(p — > 7r°e + ) ^ 5 x 10 33 yrs; 
90% C.L.) only when 32 all the low-energy uncertainties are exploited. See section I3~31 for the 
uncertainties that arise in the way of our analysis. 

29 The proton decay can be induced by the gauge-boson exchange also in SUSY-GUT models in higher- 
dimensional spacetime jS][7]. The branching ratio of various modes can be different from the standard 
one in those models. 

30 If the dimension-five decay is not the dominant process in the last class of models. 

31 We presented a numerical value of the upper bound of the proton lifetime in the abstract of this article. 
These two uncertainties are not included there, since it is impossible to make a precise estimate of them at 
this moment. The following two uncertainties, on the other hand, are included in obtaining the numerical 
value in the abstract. 

32 This statement holds as long as the non-renormalizable operator (12211 is neglected, i.e. as long asc> —1. 
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The lifetime is shortened by a factor not smaller than 0.22 if SU(5)gut-(5+5*) exists at 
low energy. The threshold corrections from these particles contribute by a factor not smaller 
than 0.47 through the change in M 2 _3, and the changes in «gut and in A R contribute by fac- 
tors not smaller than 0.66 and 0.71, respectively. Thus, those particles at 1 TeV would hardly 
be reconciled with the experimetal bound without incorporating the non-renormalizable op- 
erator 11221) . 

The lifetime is bounded from above by 

rfr- A+)< 2.1X10- (-^M 4 (10-) f°° 15GeV3 Vf— )7^ 2 y re (38) 
VF ; ~ U0 163 GeW V J\ aH J V25a GUT / \A R ) y V ; 

in the SU(5)gut x U(3)h model. Uncertainties arise 33 from c and an as in the previous 
model. The value of the QCD coupling constant is not relevant to the prediction. The SUSY 
threshold corrections changes the upper bound typically by a factor of (0.40-2.5) as M x _ 2 
changes from 10 16 20 GeV to 10 16 40 GeV. 

The lifetime is shortened by a factor not smaller than 0.56 in the presence of SU(5)gut- 
(5+5*) below the GUT scale. The decay is more enhanced as their mass is smaller. The 
enhancement factor 0.56 (when the mass is 1 TeV) consists of suppression factor 1.2, which 
comes from Mi_ 2 changed by threshold corrections of these particles, and enhancement fac- 
tors 0.66 and 0.71 respectively from ctcuT and A R . 
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A Renormalizat ion- Group Equations 

In this section, renormalization-group equations of coupling constants of the models are listed. 

SU(5)gutxU(2) h model 

33 They are not included in obtaining the numerical value in the abstract just for the same reason as in the 
previous model. 
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B J\f = 2 SUSY and Infrared-Fixed 
Renormalization-Group Flow 



Particle contents in the GUT-symmetry-breaking sector of the SU(5)gut x U(2)h model can 
be regarded as multiplets of the Af = 2 SUSY [T7], and Af = 2 SUSY is enhanced in this 
sector [Hj when the SU(5)gut gauge interaction is neglected and coupling constants satisfy 



One can see in the right panel of Fig. El that the parameter region survives in the presence of 
the 2-loop effects only when the Af = 2 SUSY is approximately preserved; the M 3 y ~ M 3 c 
line is equivalent to c*2H — o^h when the SU(5)gut gauge interaction is neglected. 

This is not a coincidence. In any gauge theory with Af = 2 SUSY, gauge coupling constants 
are renormalized only at the 1-loop level j2H]- Anomalous dimensions of hypermultiplets 
vanish |2E] at all orders in perturbative expansion, and even non-perturbatively. Therefore, 
the parameter allowed in the 1-loop analysis is still allowed when Af = 2 SUSY is preserved 
even after the 2-loop effects have also been taken into account. 

The band of the region around the (Af = 2)-SUSY line almost becomes a line as Mq 
becomes larger. For parameters above that line, the c*2H coupling constant becomes large at 
a renormalization point lower than M 3 y, while for parameters below the line, the a^H coupling 
becomes large at a renormalization point lower than M 3 c', a viable set of parameters was not 
found for large Mq even on the line in our numerical calculation. It does not mean, however, 
that the parameter does not exist at all on the (Af = 2)-SUSY line (M 3V — M 3C ) for large 
Ma-, as seen below. The (Af = 2)-SUSY relations are not only renormalizat ion-group 



9m - Am(~ A' 1H ) 



gm — -W~ Kn)- 



(52) 
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invariant but also infrared (IR)-fixed relations of the renormalization group: 





6 + A-) Kh 4 ; Qih) ,(53) 




<~ a ^ , (54) 

where the SU(5)gut interaction is still neglected. This implies that the renormalizat ion- 
group evolution to UV immediately becomes unstable 34 , even for a set of parameters that 
is slightly distant from the IR-fixed relations. The IR-fixed property (UV instability) also 
implies that the parameter region is thin only when M = M G , and not when the coupling 
constants are evaluated at M> M G - 

Thus, we can expect that the 1-loop analysis is completely reliable for a set of parameters 
exactly on the (N = 2)-SUSY line and, in particular, that a viable set of parameters does 
exist on the line even if it is not found in the numerical analysis. Therefore, the maximum 
value of M G is given at a point indicated by an arrow in the right panel of Fig. |21 At least, 
there would be no doubt that the maximum value of M G obtained in such a way provides a 
conservative upper bound of M G . 

The above argument, however, is correct only when the SU(5)gut gauge interaction is 
neglected. Therefore, let us now discuss the effects of the SU(5)gut gauge interaction. These 
break the Af = 2 SUSY in the sector. Thus, the (Af = 2)-SUSY relations in Eq. are 
no longer renormalization- group-invariant, and the renormalization-group flow is no longer 
1-loop exact. However, the SU(5)gut interaction is much weaker than the U(2) H interactions, 
and its effects are small 35 . Thus, they can be treated as small perturbations to the (A/" = 
2)-SUSY flow. In particular, the IR-fixed property of the renormalization-group equations 
(JH§j) - ()44j) is n °t changed 36 , except that the IR-fixed relations are slightly modified into 

(a^H-aifflX (tt2H-«2H)) ("lH-ain), («ih - «ih) — 0(a GJ j T ). (55) 



34 This is the reason why we believe that it would not make the parameter region wider even if we set the 

values of undetermined parameters a^H an d a m differently from those in Eq. H16(l . A deviation from the (Af 
= 2)-SUSY relation at M = Mq would immediately lead to a UV-unstable behaviour in the renormalization- 
group evolution. 

35 This can be seen from the fact that the parameter region is still almost around the (Af — 2)-SUSY line, 
i.e. M^v — M3C, in the right panel of Fig. [21 (Af = 2)-SUSY-breaking interactions are included in the figure. 

36 There is no IR-fixed relation in its strict meaning in the presence of the SU(5)gut interaction. The "IR- 
fixed relations" in Eq. Ij55(l involve o^h' ano - ^ih m ^ ne right-hand sides, and hence the "fixed relations" 
themselves change as the coupling constants flow. However, we still consider that they are almost IR-fixed 
relations, because the beta functions of the quantities in the right-hand sides (~ 0(q:gut)) are much smaller 
than those of the quantities in the left-hand sides. 
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Coupling constants flow down into the modified fixed relations and then follow the relations. 
Thus, the evolution of the coupling constants toward the UV is the most stable when the 
parameter satisfies the "IR-fixed relations" . The modified fixed relations are still almost the 
(A/" = 2)-SUSY relations, and hence the 1-loop evolution is almost correct for the parame- 
ter satisfying the relations; beta functions are different from those at one loop only by an 
order of «gut- Moreover, combinations such as (a2n - o^h ) partially absorb the SU(5)gut 
contributions in the beta functions. Therefore, the corrections to the 1-loop evolution are 
estimated conservatively from above when the SU(5)gut contribution is purely added to the 
1-loop beta functions, as we did in subsection 13.31 
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Fields 




5* X Q i ,Q i Q 6 ,Q 6 


R charges 


1 


12 



Table 1: (Mod 4)-R charges of the fields in the SU(5) GU txU(2) h model. 
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5| H(5Y 
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x a Q a i 


Q\ 




R charges 
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Table 2: (Mod 4)-R charges of the fields in the SU(5) GU txU(3) h model. 



(3,2)-« 


(1,1)° 


(1,1)° 


(l,adj.)° 


(l,adj.)° 


m.vect. 


m.vect. 


X + X ] 


m.vect. 


X + X ] 


M G = 
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M xc = 


M 3V = 


M 3C = 


V2g G uTV 


^(g 2 m + 3g 2 GVT /h)v 


V2\ m v 


\/ 2 (£2H + 3GUT> 


v / 2A 2H f 



Table 3: Summary of the particle spectrum around the GUT scale of the SU(5)gutxU(2)h 
model. The first line denotes the representation under the gauge group of the MSSM. In the 
second line, m.vect. denotes M = 1 massive vector multiplets and \ + a P arr of M = 1 
chiral and antichiral multiplets. In the last line, the mass of each multiplet are given in terms 
of gauge coupling constants and parameters in the superpotential 0. 
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M 1C = 


M 8y = 


M 8C = 


V2gv 
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^2(^ H + 2<r»/5)t; 


V2\ m v 
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Table 4: Summary of the particle spectrum around the GUT scale of the SU(5)gutxU(3)h 
model. The SU(5)gut gauge coupling constant #gut is abbreviated as g in this table. See 
the caption for Table El for the conventions in this table, replacing "superpotential ©" by 
"superpotential (J7|)". 
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Figure 1: Close-up view of the unification of the three gauge coupling constants of the MSSM. 
The fine structure constants in the DR scheme of the U(l)y, SU(2)^ and SU(3)c gauge in- 
teractions are denoted by 0^2,3, respectively. Three lines of 0J3 correspond to three different 
experimental inputs; the QCD coupling constants a^ s ^ 5 \Mz) = 0.1132 (—2a), 0.1172 (cen- 
tral value) and 0.1212 (+2<r) are used [29J. The 2-loop renormalization-group effects of the 
MSSM and the SUSY threshold corrections are taken into account. The latter corrections 
are those from the SUSY-particle spectrum determined by the mSUGRA boundary condition 
with tan/? = 10, A = GeV, (m ,mi/2) = (400 GeV, 300 GeV) and ji > (see the caption 
for Fig. H]for the convention on the sign of /i). 
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Parameter Space (1 -loop) Parameter Region (2-loop) 




Lo^ (M G /1 GeV) Log! (Mq/1 GeV) 

Figure 2: Parameter region of the SU(5)gut x U(2) h model. The parameter space of the 
model spanned by two free parameters Mq and M 3 y/M 3 Q are restricted by requiring that all 
the running coupling constants of the model remain finite while the renormalization point is 
below the heaviest particle of the model. The left panel shows the parameter region, where 
the 1-loop renormalization group is used for all the coupling constants. The right-hand sides 
of the four curves labelled "(gauge-coupling)-mass" are excluded. The region below a curve 
labelled "a^-M^c" is also excluded. Thus, the parameter space of the model is restricted to 
the shaded triangular region. The right panel shows the parameter region (shaded), where 2- 
loop effects are included in the renormalizat ion-group equations of gauge coupling constants. 
The four curves are those found in the left panel; we keep them just because they make it 
easier to compare the panel with the left one. The majority of the triangular region in the left 
panel is further excluded because of the 2-loop effects, and only a small region survives near 
the line M 3 y ~ M 3 q. The upper bound of Mq is indicated by an arrow. In the right panel, 
M 2 _3 indicates the unification point between \ja,i and l/a 3 (+2cr) (see Fig. Qfor details). 
It is easy to see that M G < 10 15 - 6 GeV ~ (l(T a4 ~ 0.40) x (M 2 _ 3 ~ 10 16 - GeV). Both two 
panels use «^ S '( 5 )(M^) = 0.1212. The effects from the non-renormalizable operator (J22j) are 
not included here. 
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Parameter region 
1 i 1 1 1 1 1 1 — 




15 15.2 15.4 15.6 15.8 16 16.2 16.4 16.6 16.8 17 
Log 10 (M^GeV) 

Figure 3: Parameter region of the SU(5)gut x U(3)h model. The parameter space of the 
model is spanned by three independent parameters: Mq, M 8 y/M$c and (Mh c Mh s )/M g . 
The figure is the \J (Mjj c Mh-) / M G = 10 03 cross section of the parameter space. We require 
that all the coupling constants remain finite under the renormalization group, while the 
renormalization point is below the heaviest particle of the model. This condition is satisfied 
in the shaded region when the 1-loop renormalization group is used. Thin curves and lines 
labelled "(gauge coupling)-mass" are lines where the corresponding gauge coupling constants 
become infinite at the corresponding mass scales. After two loop effects are included in the 
beta functions of the gauge coupling constants, the remaining allowed parameter region is 
only on the thick curve labelled 2-loop analysis. Points (A) and (B) denote the upper and the 

lower bound of the gauge-boson mass Mq, respectively, for fixed \J (Mh c M H - ) / M G = 10 a3 . 
The upper bound of Mq in the model is obtained as the maximum value Mq takes at (A) 
as \j {Mh c Mhc ) / M G changes. One also sees immediately that the lower bound at (B) is 
so low that it is of no physical importance. Mi_2 (M2-3) indicates the unification point 
between l/a\ and 1/^2 ( 1/^2 and l/a^(—2a) ) (see Fig. for more details). Note that 
(Mq at (A)_)_< 10 16 - 13 GeV ~ ((0.60 ~ 10~ a22 ) x (Mx_ 2 ^ 10 16 35 GeV)). The QCD coupling 
constant af s ' {5) {M z ) = 0.1132 is used. The effects from a non-renormalizable operator that 
corresponds to (|2*2*|) in this model are not included here. 
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Proton Life-Time (10 34 yrs.) 
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Figure 4: Contour plots of the upper bound of the proton lifetime on the mSUGRA parameter 
space. The left panel is the prediction of the SU(5)gut x U(2)h model and the right one that 
of the SU(5)gut><U(3)h model. The upper bound changes as the universal scalar mass m 
and the universal gaugino mass 7711/2 are varied (other mSUGRA parameters are fixed at 
tan/3 = 10, Aq = 0.0). The \x parameter is chosen to be positive, when the constraint 
from the branching ratio of the b — > 57 process is less severe. The upper bound of the 
lifetimes varies as (1.4-3.2) x 10 33 yrs in the SU(5)gut x U(2)h model, where the QCD coupling 
constant a^ s '^(Mz) = 0.1212 is used. The upper bound varies as (1-5) xlO 35 yrs in the 
SU(5)gutxU(3) h model, where the QCD coupling constant c^' (5) (M z ) = 0.1132 is used. 
In both panels, the effects from non-renormalizable operators are not included. The thick 
curves labelled rrih and m x are the bounds on the mSUGRA parameter space from the LEP 
II experiment in search of the lightest Higgs (m^ > 114 GeV, 95% C.L. ) |3U] and the 
lightest chargino (m x > 103.5 GeV, 95% C.L.) [2TJ. These curves are obtained by using the 
SOFTSUSY1. 7 code The excluded region changes when other codes are used; lower 

bound of mi/2 for fixed m can be higher by about 100 GeV. The code we adopt yields the 
largest pole mass of the lightest Higgs scalar among various codes available [31], and hence 
the excluded region is the smallest. 
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M 2 _ 3 /GeV (mSUGRA) 
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Figure 5: The left panels show the contour plots of the energy scale M 2 _3, where the SU(2) L 
and SU(3)c gauge coupling constants become the same. The right panels show those of the 
energy scale Mi_2, where the U(l)y and SU(2)i gauge coupling constants become the same. 
The upper panels are contour plots on the (mo, ^1/2) parameter space of mSUGRA SUSY 
breaking, the lower ones are for the (M mcss ,m S usY = ((l/24)/(47r))A) parameter space of the 
GMSB. Other parameters are fixed for both SUSY breakings; A = 0_GeV for mSUGRA 
SUSY breaking, and tan/3 = 10.0 and n > for both SUSY breakings. af s '^(M z ) = 0.1172 
is used as the QCD coupling constant in this figure. See the caption for Fig. 4 for more 
details about the region excluded by the LEP II experiments. 
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